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ABSTRACT 


Profiles of the last few Balmer lines ( H i2 -H 27^ and 
the full profiles of the HY line were computed for a grid 
of hydrogen line blanketed model atmospheres. The use of 
the last hydrogen line and the equivalent width of HY provide 
a quick means of estimating the effective temperature and 
surface gravity of B stars. Regression curves relating 
the quantum number of the last hydrogen line to the electron 
density in the model atmosphere are given. 

Furthermore, it is shown that the difference in the 
gravity determination deduced from HY equivalent widths 
using the Edmonds-Schluter-Wells formalism and the Griem 
wing formulae is, at worst, 0.18 in the common logarithm. 

Subject Headings: Atmospheres, Stellar — Early-Type Stars — 


Line Profiles 



I. INTRODUCTION 


In the course of a study of Bp stars (Klinglesmith, 1972), 
it was noted that the high hydrogen lines looked markedly 
different from those of normal main-sequence B stars. It 
was hoped that the last hydrogen line and the application 
of the Inglis-Teller (1939) formula would be an indicator 
of chemical composition peculiarities. Preliminary calculations 
indicated that this was so; although the final results, 
as described below, revealed a more complex behavior. 

II. THEORY 

The model atmospheres were taken from Klinglesmith (1971). 
They are hydrogen line (Lyman and Balmer) blanketed models 
in hydrostatic, radiative and local thermodynamic equilibrium 
covering the effective temperature range 10,000 (2,000) 

20 , 000°K and the log g range 2.5 (0.5) 4.5 for hydrogen 
mass fractions of 0.0164, 0.143, 0.667 and 1.0. 

The absorption coefficient for the high Balmer lines 

is the sum over the contributions of individual lines for 

which we use the formulae and table of TO, Y) from Griem 

(1960) with the corrections as noted by Griem (1962, 1967). 

• • 

The HY profiles were calculated using the Edmonds, Schluter 
and Wells (1967) formalism. In addition, profiles using 
the Griem wing formula were calculated for the models with 
X = 2/3. We found it convenient to rewrite the wing formula 
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as (in the notation of Griem, 1967) 


I(AX) = N I (AX)[N. /N + E + (l-E)R(N , T) ] (1) 

C -i 

where N and N. are the electron and ion densities, E is 
e i 

the ratio of the number of quasi-static perturbers to the 
number of electrons (eq. [8] of Griem, 1967) 


E = 1 + 2rr" 1/2 j\ 1/2 e" y dy, y ± = b 2 hc AX /2X 2 kT 
o 


( 2 ) 


The quasi-static line shape factor is given by 

3 /2 

I (AX) = (3TT/4) (2a 2 R OT ) (X 3 /AX 5/2 ) (b 2 -a 2 ) 
qs 


3/2 


(3) 


and the quantity 

ao-.T) - (3 2 /- 3 ) 1/2 4^72 (& f* 

(b -a ) w 

• [I + ln(AXw/AX ) ] 

3 

where a, b are the lower and upper principal quantum 
numbers, f(a, b) is given by Table 1 of Griem (1967), a Q 
is the Bohr radius, R is the Rydberg constant, 

AX = X 2 kT/chb 2 , AX_ = X 2 (e 2 N /TTmc 2 ) 1/2 and 

W r' 6 


(4) 


AX = 
a 


AX if AX < AX p 

AX if AXe[AX , AXw] 
P P 


LAX if AX > AXw 

w 


(5) 


All the other symbols have their usual meaning. The value 
of (32/tt 3 ) 1//2 is 1.032 and the quantity b 3 f(a, b)/(b 2 -a 2 ) 3//2 


is 1.0235 for HY , decreasing to 1 as b increases. 

These equations are algebraically identical to Griem 



(1960, 1962, 1967), but computationally more convenient. 

Figure 1 illustrates log 10 ){ calculated with these equations 
4 o 

for T = 10 K as a function of wavelength and electron 

pressure. The Inglis-Teller effect, which relates the upper 

quantum number of the last observable hydrogen line to 

electron density is clearly seen in this figure. The 

opacity in line center of is appreciable compared to 

its wings at log P = 0.0 and negligible at log P = 1.0 

whereas at the line opacity does not disappear until 

log P =1.5. 

& e 

III. RESULTS 

Analysis of high dispersion spectrograms showed that 
the last visible Balmer line was roughly characterized by 
its line center being depressed ~5% from the sloping "continuum” 
over the blended Balmer lines. In order to avoid correlating 
a discrete integer variable with a continuous variable (the 
electron density), we defined the quantum number of the 
last visible Balmer line as the line whose upper quantum 
number, which, interpolated , would have had exactly a 5% 
depression at line center, i.e. n(5%). We also studied 
n(20%) , but found it an insensitive parameter. 

Table 1 lists the results of the calculations. Column 
1 is the effective temperature of the model in units of 
10 3 °K, column 2 is the model log g and columns 3, 4 and 5 
are the n(5%) 's for models with X = 0.143, 2/3 and 1, 
respectively. Similarly, columns 6, 7 and 8 are the equivalent 
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widths of HY by the ESW formalism; whereas column 9 is the 
equivalent width of HY by the Griem theory for models with 
X = 2/3. Column 10 is the percentage difference between 
the ESW and Griem results (columns 7 and 9) in the sense 
(ESW - Griem) /ESW. Column 11 is the difference in log g 
from the model value if one were to fit the Griem result 
(column 9) to the "observed" ESW result (column 7). It is 
clear from columns 10 and 11 that the observational uncertainty 
is as significant as the broadening theory differences in 
determining log g. Figure 2 shows the temperature averages 
of AW(HY) and A log g. Since equivalent widths are not the 
best measure for determining the gravity, we also show the 
comparison of ESW and Griem profiles in Figure 3 for the 
16 ,000°K, log g = 4.0 model. Both the Griem and ESW 
profiles are shown as well as their difference. Similar 
plots were made of all the models referred to in Table 1; 
figure 3 is the worst case. In more than 90% of the models 
the differences between the two profiles is less than 0.025 
in the residual intensity. 

Figure 4 illustrates W(HY) versus n(5%) for the X = 2/3 
models. The lines of constant gravity and constant effective 
temperature are indicated. Similar plots can be made for 
the other values of X. The dashed lines in figure 4 are 
estimates of where the curves would go if the hydrogen line 
blanketed LTE models were valid in those regions of the T gff - 
log 6 plane. Stars for which the last hydrogen line have 
been measured (Uns51d and Struve, 1939) have been placed in 
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figure 4, where the W(HY) values are from Petrie (1953, 1965), 

Petrie and Maunsell (1950), Petrie and Moyls (1956), 

MacDonald (1953) and Williams (1936). A comparison between 

the T „„ and log g determined from this graph and the work 
eff 

of Schild, Peterson, and Oke (1971) and Heintze (1968) is 
given in Table 2. Column 1 is the star name, column 2 and. 

3 are the T and log g determined from Figure 4, column 4 

6 X X 

is Heintze 's T ff at log g = 4.0 and columns 5 and 6 are 
the T eff and log g determined by Schild, et al. One type 
of chemical composition peculiarity, namely helium under- 
abundance, is characterized by the UBV colors implying an 
earlier spectral type than the MK spectral classification. 
Table 3 lists helium poor stars (Nesterczuk, 1971) which 
are also plotted on Figure 4. It is apparent from Table 2 
and Figure 4 that strong composition anomalies can be 
detected from spectra alone. If, for example, 36 Lyn were 
really an Al star, it would lie near a Lyr and a Cma in 
Figure 4 for normal composition stars. However, HY and the 
last Balmer line indicate that, for X = 2/3, 36 Lyn should 
be about a 15,000°K, log g = 4.0 star; i.e., B5V or B6V 
which is obviously inconsistent. Figure 5 incorporates 
the visual determination of the last hydrogen line by Unsold 
and Struve (1939) with n(5%) (X = 2/3) using the Morton and 
Adams (1968) temperature scale and MK spectral types 
from the Yale Bright Star Catalog. The solid lines are lines 
of constant log g. This figure shows very clearly the surface 



gravities associated with each luminosity class. Thus 
with one spectrogram of moderate dispersion (20S/mm - 40&/mm) 
covering the region 3600$ - 4400$; the position of a star 
in figures 4 and 5 could be determined and, barring severe 
composition peculiarities, the star's effective temperature 
could be determined to within ~1000°K and its surface 
gravity to within ~0.3 in log g, which certainly is a good 
first estimate for a curve-of-growth analysis. In this 
age of fine analysis these estimates would provide a very 
reasonable starting point for a grid of models. 

In order to perform a curve-of-growth analysis the 
electron density, N , is needed rather than the surface 
gravity. It was found that the upper quantum number of the 
last hydrogen line, n(5%) correlates very well with electron 
density. Figure 6 shows the electron density at a temperature 
equal to the effective temperature for each model versus n(5%). 
The two regression equations that fit the data are given in 
the figure. Other electron densities, N g (T^) , N 6 (Tq =0.4) 
and N e (T Q = 0.1); where T R is the reversing layer temperature 
(0.8 x T eff ) and is the optical depth at the standard 
wavelength (40008) were plotted in a similar manner, however, 
N e (T eff ) produced the least variance and included the most 
models. It should be remembered that these regression 
equations are for the models. Several unpublished models 
(Klinglesmith, 1971) have been added. Six models, (20, 2.5), 
(20, 3.0), (20, 3.5), (18, 3.0), (18, 2.5) and (16, 2.5) are 
not in the regression because of their large variances. 



As long as hydrogen is the dominant continuum opacity source 
the "last visible hydrogen line" will be formed in the 
same region as the continuum since the continuum opacity 
is composed of the Paschen continuum and the merging of the 
upper levels of hydrogen. However, when hydrogen no longer 
dominates, the last hydrogen line will, in general, be 
formed in a different region from the continuum, and therefore, 
we expect little correlation between n(5%) and N e ( T e ff) 
as is the case for the low gravity, high temperature models. 
Similarly, the 10,000°K models differ from all the other 
models because H” is still a significant enough absorber 
that the last hydrogen line is formed higher in the atmosphere. 
At Tq (T ejj ) i the (10, 4.5) model has ~10 times more H~ opacity 
than the (12, 4.5) model and ~30 times more than the (10, 2.5) 
model. The latter model lies almost on the regression line 
of the hotter models. 
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IV. CONCLUSIONS 


We have shown that the last visible hydrogen line can 
be used to determine the electron density in the line 
forming region of a stellar atmosphere to within a factor 
of two. We have also shown that the last visible hydrogen 
line and the equivalent width of HY can be used to determine 
both the effective temperature and the surface gravities for 

B stars. 

Earlier investigations oif the last hydrogen line by 
van Dien (1949) and Mihalas (1966) are not really comparable 
to the present investigation. Miss van Dien used the early 
work of Pannekoek (1939) to determine gravities of stars 
observed by Yu (1926) . Tracings of Yu's observations yield 
quantum numbers consistently lower than Struve and Unsold 
by two or more except for a few higher dispersion values 
which are approximately equal to Struve and Unsold's. Her 
theoretical calculations are lower yet. Mihalas' values for 
A stars are lower by 2 for main sequence stars and by 4 for 
giants than those of the present work. His calculations, 
however, are based upon models without Lyman line blanketing 
and therefore are not expected to be comparable. 

It should be pointed out that this analysis was done 
with LTE, hydrogen line blanketed models. The effects 
of non-LTE and/or metal line blanketing have not been 
considered and remain to be done. 
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TABLE 1 

Last Hydrogen Line and HY Equivalent 


o n 0%) 

Te/10 J log g X=0 . 143 2/3 1 X=0.143 

10 2.5 22.666 24.142 24.697 8.52 

3.0 21.013 22.370 22.784 11.30 

3.5 19.642 20.745 21.130 14.60 

4.0 18.289 19.216 19.641 18.50 

4.5 16.863 17.860 18.245 22.70 

12 2.5 21.629 23.163 23.830 6.39 

3.0 20.106 21.526 22.024 8.39 

3.5 18.689 19.947 20.412 10.80 

4.0 17.366 18.506 18.877 13.60 

4.5 16.144 17.123 17.595 17.10 

14 2.5 20.887 22.586 22.968 5.06 

3.0 19.448 20.802 21.274 6.88 

3.5 18.088 19.272 19.771 8.36 

4.0 16.805 17.908 18.371 11.50 

^ 4.5 15.656 16.711 16.987 14.40 

© 16 2.5 20.502 21.823 22.421 3.74 

3.0 18.810 20.195 20.627 5.41 

3.5 17.584 18.737 19.162 7.34 

4.0 16.311 17.175 17.827 9.68 

4.5 15.167 16.236 16.580 12.40 

18 2.5 19.226 21.194 21.880 2.46 

3.0 18.166 19.698 20.193 3.89 

3.5 16.936 18.347 18.706 5.63 

4.0 15.897 16.902 17.446 7.76 

4.5 14.801 15.827 16.170 10.30 

20 2.5 18.458 20.799 21.540 1.77 

3.0 16.849 19.122 19.712 2.68 

3.5 16.218 17.838 18.284 4.15 

4.0 15.199 16.648 16.916 5.89 

4.5 14.368 15.482 15.823 8.10 


from Model Atmospheres- 


w esw^ hy ^ 


W q (HY, 


2/3 

1 


AW 

A log c 

6.28 

5.66 

6.79 

-0.081 

0.110 

8.60 

7.83 

9.30 

-0.081 

0.125 

11.40 

10.50 

12.30 

-0.079 

0.128 

14.90 

13.80 

15.90 

-0.067 

0.125 

18.90 

17.70 

20.00 

-0.058 

0.137 

4.66 

4.15 

5.05 

-0.084 

0.113 

6.39 

5.83 

6.94 

-0.086 

0.131 

8.49 

7.80 

9.19 

-0.082 

0.152 

10.80 

10.00 

11.60 

-0.074 

0.138 

13.7 0 

12.80 

14.70 

-0.073 

0.172 

3.68 

3.24 

3.99 

-0.084 

0.102 

5.20 

4.74 

5.67 

-0.090 

0.131 

7.00 

6.43 

7.60 

-0.086 

0.146 

9.05 

8.40 

9.78 

-0.081 

0.148 

11.50 

10.70 

12.40 

-0.078 

0.183 

2.90 

2.57 

3.14 

-0.083 

0.083 

4.34 

3.95 

4.73 

-0.090 

0.123 

5.93 

5.50 

6.46 

-0.089 

0.151 

7.68 

7.23 

8.32 

-0.083 

0.139 

9.99 

9.30 

10.70 

-0.071 

0.154 

2.20 

1.99 

2.37 

-0.077 

0.061 

3.60 

3.34 

3.93 

-0.092 

0.109 

5.11 

4.47 

5.57 

-0.090 

0.139 

6.76 

6.37 

7.33 

-0.084 

0.139 

8.81 

8.26 

9.48 

-0.076 

0.163 

1.70 

1.45 

1.82 

-0.076 

0.045 

3. 02 

2.81 

3.29 

-0.089 

0.102 

4.35 

4.17 

4.75 

-0.092 

0.128 

5.91 

5.64 

6.41 

-0.085 

0 . 134 

7.77 

7.40 

8.37 

-0.077 

0.161 




TABLE 2 


Comparison of Effective Temperatures and 
Surface Gravities 


Star 

T e /10 3 

log g 

T /10 3 
e 

Heintze 

T /10 3 
e 

Oke , Schild, 

log g 
Peterson 

V Ori 

27.2°K 

4.5 

29 . 1°K 



T Sco 

28.4 

4.41 

25.8 



cp Ori 

25.6 

4.18 




P Cep 

23.8 

3.85 


24.2°K 

4.0 

e Per 

23.8 

3.85 

25.4 

28.6 

4.0 

Y Peg 

20.8 

3.77 

20.4 

21.7 

4.0 

P CMa 

26.0 

3.3 




102 Her 

24.0 

3.05 




e CMa 

26.0 

2.95 


28.6 

4.0 

p Leo 

26.7 

2.83 




h Cas 

22.5 

2.54 




9 Cep 

19.3 

2.4 




C Ca s 

18.0 

3.85 




t Her 

16.2 

3.9 

14.7 



i Her 

17.8 

3.6 

16.8 



e Cas 

17.6 

3.6 


15.9 

4.0 

n T 

i n e 

O OK 




•l Juy jl 

jl i • 





C Dra 

18.0 

2.85 

13.2 



X 2 Ori 

17.0 

2.1 




X Aur 

16.0 

2.2 




P Ori 

14.2 

2.3 




55 Cyg 

15.0 

2.2 




4 Lac 

12.0 

<2.0 




e Del 

14.2 

3.6 




67 Oph 

16.0 

2.6 




a Dra 

10.1 

3.75 




a Cma 

9.8 

4.0 


10.2 

4.35 

a Lyr 

9.7 

3.9 


9.65 

4.05 
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TABLE 3 


LAST HYDROGEN LINE AND HY EQUIVALENT WIDTHS 
OBSERVED IN HELIUM- WEAK STARS 


Star 

Spectral 

Type 

UBV 

Spectral 

Type 

W<HY) 

n (5%) 

Te/10 3 

log g 

He/H 

HD21699 

B8 III 

B5 IV 

7.4 

18.0 

15.8 

3.98 

0.008 

20 Eri 

AO Si 

B6 V 

8.4 

17.7 

14.7 

4.25 

0.003 

22 Eri 

B8 

B5 V 

7.4 

18.6 

14.8 

3.95 

0.008 

20 Tau 

B7 III 

B8 III 

7.0 

20.2 

13.0 

3.45 

0.08 

12 CMa 

B7 III 

B4 V 

7.4 

17.2 

17.7 

4.43 

0.016 

36 Lyn 

A1 

B7 IV 

7.8 

19.3 

13.4 

3.82 

0.01 
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FIGURE TITLES 


1. Log^Q* (High Balmer lines only) versus log-^Q and X (8) 
for T = 10 4 °K. 

2. Temperature averages of Aw(HY) and -Alog g from Table 1. 

3. Profiles of HY from the 16,000°K, log g = 4.0 model showing 
ESW and Griem as upper and lower curves. The difference 

is shown at the bottom. Note that the AR scale is expanded. 

4. n(5%) versus W(HY) with stellar observations. Different 
luminosity classes are marked by separate symbols. Solid 
dots are the helium-weak stars. 

5. n(5%) versus spectral type with the Struve and Unsold data 
using MK spectral types and the Morton-Adams temperature 
scale. Different luminosity classes are marked by 
separate symbols. 

6. Log N q (T e ^ f ) versus log n(5%) from the models. 
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